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We have investigated correlated electronic structures and the phase diagram of electron-doped 
hydrocarbon molecular solids, based on the dynamical mean- field theory. We have found that the 
ground state of hydrocarbon-based superconductors such as electron-doped picene and coronene is 
a multi-band Fermi liquid, while that of non-superconducting electron-doped pentacene is a single- 
band Fermi liquid in the proximity of the metal-insulator transition. The size of the molecular 
orbital energy level splitting plays a key role in producing the superconductivity of electron-doped 
hydrocarbon solids. The multi-band nature of hydrocarbon solids would boost the superconductivity 
through the enhanced density of states at the Fermi level. 

PACS numbers: 71.27.+a, 74.70.Kn, 74.70.Wz, 74.20.Pq 
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Since the discovery of high Tc superconductors, the 
role of the electronic correlation in the superconductiv- 
ity has been a subject of intensive investigation. While 
conventional BCS superconductors such as Nb and MgB2 
have good metallic nature of Fermi liquid, [1, 2] unconven- 
tional superconductors such as doped cuprate and iron 
pnictide show bad metallic behavior due to their strong 
electronic correlation. [3, 4] The correlation issue exists in 
carbon-based 7r-electron superconductors too. Ca doped 
graphite (CaCe) shows conventional superconductivity 
with weak electronic correlation, [5] whereas CsaCeo ex- 
hibits unconventional superconductivity with strong elec- 
tronic correlation. [6| 

New TT-electron superconductors have been recently 
discovered in polycyclic aromatic hydrocarbon (PAH)- 
based molecular solids: Kapicene (Tc=18 K), Kacoronene 
(rc=15 K), Kaphenanthrene ffc=5 K), and K3l,2;8,9- 
dibenzopentacene (Tc=33 K).J7H1Q| These superconduc- 
tors were also reported to have strong correlation. |lll- 
Il4j On the other hand, a similar PAH-based molecular 
solid, K-doped pentacene, does not have superconduc- 
tivity, but exhibits only the metal-insulator transition 
(MIT) behavior. |l5| Note that both picene and pentacene 
are composed of five benzene rings with slightly different 
arrangements, as shown in Fig. [H The different ground 
states in K-doped picene and K-doped pentacene sug- 
gest that the correlation effects come into play distinctly 
between superconducting and non-superconducting sys- 
tems. 

In this Letter, in order to resolve the issue of cor- 
relation effects in hydrocarbon-based superconductors, 
we have investigated their electronic structures sys- 
tematically, employing the dynamical mean-field theory 
(DMFT). Based on the ground state electronic struc- 
tures, we have constructed the phase diagram of hy- 
drocarbon molecular solids as functions of doping and 



relevant energy parameters including the Coulomb cor- 
relation, the Hund coupling, and the molecular-orbital 
(MO) energy level splitting. Our studies reveal that 
hydrocarbon-based superconductors belong to multi- 
band Fermi hquid system, while non-superconducting K- 
doped pentacene belongs to single-band system in the 
proximity of the MIT. Further, we have shown that 
the energy level splitting between LUMO+1 and LUMO 
(lowest unoccupied MO) plays a key role in the super- 
conductivity of electron-doped hydrocarbon molecular 
solids. 

Both picene and pentacene solids have layered crystal 
structures, as shown in Fig. [TJa). Picene molecule has 
an arm-chair type structure (Fig. [T]^b)), while pentacene 
molecule has a linear arrangement of benzene rings 
(Fig. [TJc)). The difference in molecular structures pro- 
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FIG. 1: (Color online) (a) Crystal structure of picene and 
pentacene solids. Blue ellipses represent picene or pentacene 
molecular units shown in (b) and (c). ti, t2, ts are hopping 
amplitudes along each direction, (b) Molecular structure of 
picene. (c) Molecular structure of pentacene. (d) MO en- 
ergy levels and occupation of electrons for picene^" solid, (e) 
The same for pentacene^" solid. Note that the energy level 
splitting (A) between LUMO+1 and LUMO for picene sohd 
(A =0.040 eV) is much smaller than that of pentacene solid 
(A =1.260 eV). 



TABLE I: Hopping amplitudes ti,rn,R (in meV) determined from the ab initio band structures of picene and pentacene solids, 
where / and m represent the MOs, and R represents the hopping direction in the lattice. Subscript 1,2,3 represent the hopping 
directions as shown in Fig[T] and c represents the direction normal to ab plane. 
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duces the different electronic structures in three electron- 
doped picene (picene^") and pentacene (pentacene^"). 
As shown in Fig. [TJd) and (e), the energy level split- 
ting A between LUMO+1 and LUMO of picene solid 
(A = 0.04 eV) is much smaller than that of pentacene 
solid (A=1.26 eV). Because the bandwidth of each or- 
bital in both molecular solids is ~0.25 eV, picene^" has 
the occupation of three electrons on nearly two-fold de- 
generate orbitalsjlli llSl Il6| whereas pentacene^" has 
one electron on the single LUMO+1 orbital that is far 
separated from the lower LUMO. 

The correlation effects in electron-doped hydrocarbon 
solids are dealt with by the following two-band Hubbard 
model Hamiltonian, 

H = Hq -\- Hj^ 

= 2_^ ti^rn,IlCiji^^Cm,0,a^Hi, (1) 

where Hq and Hj are non-interacting and interacting 
Hamiltonians of doped electrons, respectively. Here 
ti,m,ii corresponds to hopping from (0,7ti) to (R,/), where 
and R represent sites, m and / represent the MOs. We 
have determined t^,rn,R, using the downfolding scheme 
of Kohn-Sham orbit als in the maximally localized Wan- 
nier function (MLWF) basis. 17H19| Kohn-Sham orbitals 
were obtained in the generalized gradient approxima- 
tion (GGA), by employing the full-potential augmented 
plane wave (FLAPW) band method [20[ implemented in 
WIEN2k package. [2 Ij For electron doped systems, we 
have utilized the rigid band approximation due to the 
absence of experimental crystal structures. The obtained 
hopping parameters ti^rn.Yi are provided in Table [H 
Hi is given approximately by 



m,R m>^,R,(j 

+ {U' -J) ^ n^,R,an^,R,a, (2) 

7n>^,R,cr 

where [/, U' and J are intra-, inter-orbital Coulomb cor- 
relation and Hund interaction parameters, respectively. 
We considered here rotationally symmetric interaction, 
so that Ulumo+1 ^ Ulumo and U' ^ U — 2 J. [13] 
We have solved the above Hamiltonian by carrying out 
the DMFT calculation. [22| We used the continuous time 
quantum Monte-Carlo (CTQMC) method as an impurity 
solver. [23, ,M] We set temperature at T=77 K (=6.67 



meV), which is low enough to observe the MIT in the 
phase diagram of hydrocarbon solids. [25] 

Figure [2] shows the imaginary part of the obtained 
Green's function G{iuJn) and the self-energy S(zco'n) for 
LUMO+1 and LUMO of picene^- solid. For U = 0.45 
eV, lmG{iuJny^ of both LUMO+1 and LUMO have finite 
values at cjn = {un- Matsubara frequency), implying 
that the densities of states (DOSs) at the Fermi level 
{Ep) are finite for both orbitals. Also ImI](ico'n)'s for 
both MOs converge to zero in the zero frequency regime, 
signifying the Fermi liquid nature. Note that the mag- 
nitude of lm.Ti{iujn) is significantly larger for LUMO+1 
than for LUMO, which reflects that the electronic cor- 
relation is stronger for LUMO+1. In fact, U = 0.45 eV 
corresponds to the value obtained from the first princi- 
ples calculation for picene^" solid. [13] So these results 
indicate that Kapicene solid in its normal state has two- 
band Fermi liquid nature with clear orbital-selective band 
renormalization. 

For [7=0.65 eV, ImG(ia;n)'s near Un = appear to be 
finite for both LUMO+1 and LUMO, but that for LUMO 
has negative slope. From the derivative oflinG{iuJn) near 
cjn = one can distinguish between metallic and insu- 
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FIG. 2: (Color online) Imaginary Green's functions G{iuJn) 
and self energies T^{iuJn) of LUMO+1 and LUMO for picene^" 
solid (J=0.05 eV and A=0.04 eV). Red (dotted), green 
(dashed), and blue (solid) lines are for [7=0.45 eV, [7=0.65 
eV, and [7=1.00 eV, respectively. 
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FIG. 3: (Color online) (a) Ground states of picene solid 
(J=0.05 eV and A=0.04 eV) with variations of doping 
(N=l,2, and 3) and Coulomb interaction (U). (b) Ground 
states of pentacene solid (J=0.05 eV and A=1.26 eV). IFL, 
2FL, BI, and MI stand for one-band Fermi liquid, two-band 
Fermi liquid, band insulator, and Mott insulator, respectively. 
Dashed line for picene^" corresponds to the U value obtained 
from the first-principle calculation. [13.] 



lating phases as shown in previous literatures. |26i |27|. 
We consider that positive (negative) derivative illustrates 
the metallic (insulating) electronic structure. Accord- 
ingly, only the LUMO+1 band has nonzero DOS at Ep. 
The behaviors of lm.G{iuJn) and lmE{iuJn) of LUMO+1 
for t/=0.65 eV reflect that the system belongs to nar- 
row single-band Fermi liquid state. For [7=1.0 eV, both 
ImG(zco'n)'s are vanishing at ujn = 0, implying that DOSs 
at Ef are zero for both LUMO+1 and LUMO. lmJ:{iuJn) 
at cjn = diverges for LUMO+1, indicating that the sys- 
tem has the Mott insulating state with hole-orbital dis- 
proportionation nature at LUMO+1. Namely, the hole 
orbital nature changes from the mixture of LUMO and 
LUMO+1 to the single LUMO+1 (see the inset of Fig.©. 
Figure [3] shows the ground states of picene and pen- 
tacene solids depending on doping N and the U value. 
We have obtained the ground states by following the 
steps in Fig. [21 As mentioned above, picene^" solid, 
which has U = 0.45 eV, belongs to the two-band 
Fermi liquid state. |28| This feature reveals that the su- 
perconductivity in three electron-doped picene, such as 
Kapicene and Cai.spicene, emerges from the multi-band 
Fermi liquid state. If we increase U further for picene^", 
single-band (LUMO+1) Fermi liquid state is realized at 
[7=0.625 eV, and the Mott insulating state at /7=0.90 
eV. The occurrence of the MIT for large U is consis- 
tent with the previous finding. |14| Similarly to picene^", 
picene^" solid also shows the phase transition upon in- 
creasing /7, from two-band to one-band (LUMO) Fermi 
liquid state, and then to the Mott insulating state. In 
the case of picene^", due to the finite Hund coupling J 
(=0.05 eV), the system for small U remains as two-band 
Fermi liquid state with high-spin configuration. With in- 
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FIG. 4: (Color online) (a) Quasi-particle residue Z's of 
LUMO+1 and LUMO vs. U for picene^" soUd with varia- 
tion of A (J=50 meV). (b) The same for picene^" solid with 
variation of J (A=40 meV). Red arrow denotes the difference 
of Z's between LUMO+1 and LUMO. 



creasing [7, the phase transition occurs directly from two- 
band Fermi liquid to Mott insulating state. If one takes 
into account the non-rigidity of band structures arising 
from the hybridization with cations in K or Ca doped 
picene, the electronic correlation effect would be further 
reduced with respect to the above rigid band case. [19) 
Then Kapicene and Cai.spicene would have more stable 
two-band Fermi liquid nature in their normal states. 

Noteworthy in Fig. ^h) is that the ground state of 
pentacene^" is very different from that of picene^" . Due 
to much larger A value in pentacene, both pentacene^" 
and pentacene^" solids exhibit the transition from single- 
band Fermi liquid to Mott insulator, like a single-band 
half- filled system. On the other hand, pentacene^" has 
the low-spin state due to large A, and so only the band 
insulating state is realized. The findings in Fig. [3] mani- 
fest that different A's of picene and pentacene play key 
roles in producing the different phases. Also high-spin 
and low-spin configurations of picene^" and pentacene^" 
suggest that the Hund couphng J affects the phase dia- 
gram. 

Figure |4] presents the quasi-particle residue Z for 
picene^" solid with the variation of A and J. The quasi- 
particle residue Z is obtained from Z = (1 — ^...!^° )""*", 
where wq is the lowest Matsubara frequency. [26] With 
increasing A, the following features are observed in 
Fig. m^a): (i) the larger difference of Z's between 
LUMO+1 and LUMO, (ii) the wider single-band 
(LUMO+1) Fermi liquid regime, and (iii) the reduced 
critical U for the MIT. The above observations indi- 
cate that just a small enhancement of A would induce 
the single-band Fermi liquid state in picene^" sohd, 
and thereby suppress the superconductivity. This sensi- 
tive dependence of electronic structure of picene^" upon 
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FIG. 5: (Color online) Phase diagram of hydrocarbon^" solids 
with respect to Coulomb correlation U and MO level splitting 
A. Positions of picene^", coronene^", and pentacene^" solids 
are marked in the phase diagram. Inset diagram represents 
the schematic DOS of each phase. In the MI regime, the 
hole-orbital disproportionation nature is realized. 



variation of A explains different experimental electronic 
structures of K epicene solid depending on the prepa- 
ration condition. jl4 1291431] As shown in Fig. IH^b), the 
larger J also yields the larger difference of Z's between 
LUMO+1 and LUMO. But the single-band Fermi liq- 
uid regime and the critical U value for the MIT do not 
change much. Lager difference of Z's between two MOs 
implies that the inter-orbital fluctuation becomes sup- 
pressed. Therefore Fig. |4] indicates that the larger A 
induces the orbital disproportionation, and the larger J 
suppresses the inter-orbital fluctuation. |32| 

Figure [5] presents the overall phase diagram of 
hydrocarbon^" solids as functions of U and A. The 
ground state of picene^" solid is two-band Fermi liquid 
state. In contrast, the ground state of pentacene^" solid 
is single-band Fermi liquid or Mott insulator depending 
on U value. [33] The ground states of Kacoronene and 
Ksphenanthrene in their normal states are also two-band 
Fermi liquid states, with three doped electrons occupy- 
ing nearly degenerate LUMO and LUMO+1. On the 
other hand, l,2;8,9-dibenzopentacene has large energy 
splitting between LUMO and LUMO+1, but nearly de- 
generate LUMO+1 and LUMO+2.[33 Hence K3l,2;8,9- 
dibenzopentacene has one electron on the nearly degen- 
erate LUMO+1 and LUMO+2. This situation is similar 
to that of picene^" in Fig. [3fa). Therefore the phase 
diagram of Fig. [5] demonstrates that the superconduc- 
tivity in hydrocarbon-based superconductors commonly 
emerge from the two-band Fermi liquid state. The under- 
lying superconducting mechanism in hydrocarbon-based 
superconductors needs further investigation, but it is ev- 
ident that the superconductivity would be boosted up by 
the enhanced DOS at Ej? due to their two-band Fermi 
liquid nature. [35j 



In conclusion, we have investigated the electronic 
structures of electron-doped hydrocarbon solids based on 
the DMFT calculations, and constructed the phase dia- 
gram with respect to Coulomb correlation [/, doped elec- 
trons A/", Hund coupling J, and MO energy level splitting 
A. We have shown that the superconductivity in hy- 
drocarbon solids commonly emerges from the two-band 
Fermi liquid state. This is in contrast to the case of 
non-superconducting pentacene, which has the effective 
single-band Fermi liquid state in the proximity of the 
MIT. The size of MO energy level splitting plays an im- 
portant role in determining the ground states of hydro- 
carbon solids. Our results demonstrate that multi-band 
nature in hydrocarbon solids is essential to boost the su- 
perconductivity through the enhanced DOS at Ep- It 
is thus suggested that higher Tc superconductors need 
to be searched for in molecular solids with multi-band 
nature at Ejr. 
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